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Abstract 

We discuss a model for the lepton sector based on the seesaw mechanism and on 
a -D4 family symmetry. The model predicts the mixing angle #13 to vanish. The solar 
mixing angle #12 is free — it will in general be large if one does not invoke finetuning. 
The model has an enlarged scalar sector with three Higgs doublets, together with 
two real scalar gauge singlets %i (* = 1)2) which have vacuum expectation values 
{Xi)o a t the seesaw scale. The atmospheric mixing angle #23 is given by tan #23 = 
(X2)o/(Xi)n) arm ^ is maximal if the Lagrangian is L>4-invariant; but D4 may be 
broken softly, by a term of dimension two in the scalar potential, and then (x2)o 
becomes different from (xi)o- Thus, the strength of the soft D4 breaking controls the 
deviation of #23 from 7r/4. The model predicts a normal neutrino mass spectrum 
(m.3 > 7712 > fn{) and allows successful leptogenesis if m\ ~ 4 x 10~ 3 eV; these 
properties of the model are independent of the presence and strength of the soft D4 
breaking. 
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1 Introduction 



The idea of neutrino oscillations pQ through neutrino masses and lepton mixing j2] has 
proved successful and has led to continuous theoretical and experimental efforts and 
progress — for a review see, for instance, [3]. Still, no favoured theory for explaining the 
observed patterns of neutrino masses and lepton mixing has emerged yet; in recent years, 
myriads of textures and models for the lepton mass matrices have been proposed — for a 
review see, for instance, [3] . Only the smallness of neutrino masses has found a favourite 
explanation in the seesaw mechanism j^]. In this paper we employ that mechanism, and 
dismiss textures, trying instead to trace some features of the lepton sector to family sym- 
metries. In particular, one may try and elucidate whether there is a connection between 
the smallness of the mixing angle # 13 and the maximality of the atmospheric mixing angle 
#23- It turns out that there is no such connection: 

• In the model of jE], which is based on the non-Abelian symmetry group 0(2) j7j, 
and in the model of jH], which is based on the discrete symmetry Z) 4 , one finds 
#13 = and #23 = tt/4; the model of S., which is based on the permutation group 
S3, displays only a minimal deviation from those predictions. 

• The authors of JU], on the basis of a symmetry A4, found #23 — 7r /4 but their # 13 
is non-zero and can even be rather large. 

• With the non-standard CP transformation used in ^1] atmospheric mixing is max- 
imal but #i 3 is free. 

• Conversely, it was shown in ^21 that one can construct models based on Abelian 
groups, in particular on Z4, which display #13 = but have a free atmospheric 
mixing angle. 

The listing above demonstrates that, without specifying the family symmetry group and 
other details of the model, no general statement can be made as regards a possible relation 
between #13 and #23- 

In this paper we want to discuss a modification of the D4 model of 8\. In this 
modification one breaks the -D4 symmetry softly in such a way that #13 = is preserved 
but the atmospheric mixing angle becomes free and controlled by the strength of the 
soft symmetry-breaking term. The model discussed in this paper has the same mixing 
features as the models in [TJ] , but a completely different structure due to the non-Abelian 
symmetry group. Both the original and the softly broken D 4 models belong to a class of 
models where the Yukawa-coupling matrices are diagonal, and lepton mixing originates 
exclusively in the Majorana mass matrix M# of the heavy neutrinos; the suppression 
of neutral flavour- changing interactions is a general feature of that class of models jT3]. 
We furthermore show that, in the softly broken D 4 model, successful leptogenesis can be 
traced back to the original D4 model. 

This paper is organized as follows. In Section El we review the original D4 model. The 
softly broken D4 model is introduced in Section El where we also investigate its predic- 
tions for lepton mixing. The neutrino mass spectrum, the effective mass in neutrinoless 
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(3(3 decay, and leptogenesis are studied in Section EJ In Section we perform a qualita- 
tive discussion of the renormalization-group corrections to the prediction # 13 = 0. Our 
conclusions are found in Section |3 



2 The original model 

The -D4 model [S] has three lepton families, including three right-handed neutrinos which 
enable the seesaw mechanism We generically denote e, /1, and r by a; we have then 
three left-handed lepton doublets D aL , three right-handed charged- lepton singlets a R , 
and three right-handed neutrino singlets v a R. The scalar sector consists of three Higgs 
doublets 4>i, 2 , and <p 3 , together with two real gauge singlets Xi an d X2- 
The model has three symmetries of the Z2 type: 



4 T) 

r^(aux) 



D tL , t r , v rRy X2 change sign; 

D^ L <-+ D tL , /Ir <-+ Tr, V^r <-> V t r, Xl 

VeR-, v^r, v tR , 0i, e R change sign. 



X2, 03 



The symmetry zjj flips the sign of all multiplets of the r family, while l^ 1 ' exchanges 
the multiplets of the /1 and r families. The auxiliary symmetry Z 2 aux ' ) ensures that 02 and 

(tr) 

03 do not have Yukawa couplings to the v aR , so that the spontaneous breaking of 7L\ , 
which is necessary for m M 7^ m r , does not have consequences in the neutrino Dirac mass 
matrix Mo- 

As discussed in jS], Z 2 T ^ an d Z^ do not commute, and together they generate a 
non-Abelian group D4 with eight elements. This group has five inequivalent irreducible 
representations (irreps): one two-dimensional irrep and four one-dimensional irreps. It 
is clear that (D^, D tL ), (fi R ,r R ), (u^r, v t r) , and (xi,X2) transform as doublets of _D 4 , 
whereas 0i and 02 transform according to the trivial one-dimensional irrep, and 03 ac- 
cording to a non-trivial one-dimensional irrep. 

The above multiplets and symmetries determine the Yukawa Lagrangian 

- yiD eL u eR + y 2 [p^Vfjjt + D rL v TR j 0i 

-y 3 D eL e R (j)i - 2/4 (p^L^n + DtlTr) 02 - y 5 [D^m - D tL trJ 3 



(tr) 



C 



Y 



+ 2 Vx v eR C 1 (*V*Xl + VtRXz) + H.c, 



where 0i = iT%(j)\; they also determine the Majorana mass terms 
\ [M*v T eR C~ l v eR + M'* (^ r C-\r + v t tR C~ x 

'0) = vJV2 (j = 



M 



H.c. 



(2) 



(3) 



The vacuum expectation values (VEVs) ^0 
lepton masses as 

\f2m e = |j/ 3 t>i|, 

v^m^ = 12/4^2 + 2/5^1 

v / 2m T = 1 2/4^2 - 2/5^3 1 



1, 2, 3) fix the charged- 



(4) 
(5) 
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The neutrino Dirac mass matrix is given by 

M D = diag (a, b, b) , (6) 

with \/2a = y*Vi and \p2b = y^vi- The charged- lept on mass matrix and Mn are simulta- 
neously diagonal. Therefore, lepton mixing must result exclusively from a non-diagonal 
M R , and this arises by virtue of the non-zero VEVs of xi,2- Let us write those VEVs as 

(0| X i|0) = Wcos 7 , 

(0 |Xa| 0> = Wsinj, 1 j 

with W real and positive. Then the Majorana mass matrix for the right-handed neutrino 
singlets is 

/ M M x cos7 M x sin7 
Mr (7) = M x cos7 M' 
\ M x sin7 M' 

where M x = y*W. The effective Majorana mass matrix for the light neutrinos is given by 

M v (rf) = -M T D M R { 1 y l M D , (9) 

where Me, = is in equation (|5jl. Since the charged-lepton mass matrix is diagonal, 
the lepton mixing matrix U is found simply through the diagonalization procedure 

U T M U (7) U = diag (mi, m 2 , m 3 ) , (10) 

with real and non- negative masses rrij (j — 1, 2, 3) for the light neutrinos. 
In the .D 4 -invariant scalar potential V sym , the x-dependent terms are jH] 

^sym = ■ • • + (Xl + xl) ["A* + P10101 + P2 (0202 + 4<fo)] + A (*? + ' 

+ (xl ~ Xl) (v4fo + ri*4fa) + A' (xl -xl)\ (11) 

where /i, px, p 2) A, 77, and A' are c-numbers. Then, 7 is determined by the minimization 
of 

/ sym (7) = \'W 4 cos 2 2 7 + Re (77^3) W 2 cos 2 7 . (12) 
Provided A' > and |Re (riv^y^)] < 2AW 2 , the minimum of / is at 

Re {r]v^v 3 ) , * 

COS27 = -^AW^- (13) 

According to the seesaw mechanism j^j, we assume 

|M|,|M'|>t>, (14) 



where v = y |t>i| 2 + |f2| 2 + | ^3 1 2 — 246 GeV. We furthermore assume W ~ \M\ , \M'\, so 
that the off-diagonal matrix elements of Mr (7) are not much smaller than the diagonal 
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matrix elements, lest the solar mixing angle becomes very small. Under these assumptions, 
|cos27| is negligibly small and, for all practical purposes, 



<0| X i|0> = <o| X2 |o> 



w_ 

71' 



(15) 



or 7 = 7r/4. Corrections to this value of 7 are of order v 2 /W 2 . 
One can easily check that 



SM R {n/4)S = M R {n/4), 



(16) 




(17) 



where 

/ 1 

S = 

\o 

Then, due to equations © and ©, also SA4 U (7r/4) S = A4 u (ir/4). From this the 
predictions of the D4 model follow 8J: maximal atmospheric mixing angle, i.e. 623 = 7r/4; 
free, in general large, solar mixing angle 9±2] and 613 = 0. One also finds — see Section^]— 
that, because the (/x, r) matrix element of M R (7) vanishes, the inverted neutrino mass 
spectrum is excluded. 



3 The softly broken D4 model 

In the original -D4 model of the previous section, all three symmetries of equation (0) 
are spontaneously broken, and none of them is broken in the Lagrangian itself. In this 
section we introduce a soft dimension- two breaking of Z 2 . It is easily found that there 
are two terms performing that breaking: one of them is </>203 plus its Hermitian conjugate, 
which is largely irrelevant; the other one is 

Koft = //soft (x? - xl) > ( 18 ) 

which is to be added to the scalar potential. Now the angle 7 of equation (JJJ) is determined 
by the minimization of 

/(7) = /sym(7) +/i S oft^ 2 COs2 7 , (19) 

leading to 

// soft + Re (77^3) 
cos 2 7 = • (20) 

We assume |/J so ft| to be of the same order of magnitude as W, i.e. we assume it to be of 
the seesaw scale; this results in a non-negligible deviation of 7 from 7r/4. As discussed in 
the previous section, such a deviation is in general present, cf. equation (|T3|) . yet without 
the soft breaking in equation (|TH|) it is very small because v <C W. 1 

1 We might as well assume W to be of the Fermi scale, and then /i so f t would be unnecessary. But this 
would imply the off-diagonal elements of Mr being much smaller than the diagonal ones, and this would 
prevent the fitting of the observed large solar mixing angle. 
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Defining 

/ 1 \ 

5(7)= cos2 7 sin2 7 = S (7)^ , 

\ sin 27 — cos 27 J 

which generalizes the matrix S in equation (|17jl. it is easy to check that 



(21) 



(22) 



S( 1 )M R ( 1 )S( 1 ) = M R ( 1 ) 5( T )M i/ (7)S , (t)=^(t) 
The orthogonal matrix S 1 (7) has a unique eigenvalue —1, corresponding to the eigenvector 

/ 



u 3 



— sin 7 

y cos 7 



(23) 



Equation (J22|l then implies that u 3 is also an eigenvector of M. v (7). Since u 3 is real and 
its first entry is zero, this eigenvector — possibly multiplied by a phase — must constitute 
the third column of the lepton mixing matrix U. Explicit calculation of Ai u (7) using 



1 



M' (MM' - Ml 



( 



M' 2 -M'M X cos 7 

-M'M X cos 7 MM' - Ml sin 2 7 



-M'M X sin 7 
M 2 sin 7 cos 7 



(24) 



MM' - M 2 cos 2 7 



\ — M'M X sin 7 M 2 sin 7 cos 7 
shows that the eigenvalue is —b 2 /M', hence the mass of the third light neutrino is 

b 2 



Thus, 



is exact at the tree level, while 



m 3 



U, 



M' 



c3 



sin 2 26. 



atm 



4 \U, 



(i3 I 



It/, 



(i3 I 



sin 27, 



(25) 



(26) 



(27) 



i-e- ^atm — 7- The atmospheric mixing angle is equal to the angle 7 from the VEVs of 
Xi,2 in equations (J7J). 

We thus have a model where the atmospheric mixing angle is arbitrary, since its 
deviation from 7r/4 is controlled by the strength /i soft of the soft-breaking term V^ f t in 
equation (fTHj) ; on the other hand — at least at the tree level — the angle #13 vanishes. The 
solar mixing angle is completely arbitrary, but in general it will be large. 

We have broken the /x-r interchange symmetry softly in the scalar potential by 
means of the term in equation (|18J) . That soft-breaking term is of dimension two. One 
might, in addition, also break 7L 2 through terms of dimension three, viz. by assuming 
(Mr) to be different from (Mr) tt . A general mixing matrix U would then follow. It is, 
however, consistent to avoid soft-breaking terms of dimension three while allowing those 
of dimension two. 



6 



4 Neutrino masses and leptogenesis 

The diagonalization of the mass matrices in equations (JHJ) and Q can be related to the 
diagonalization of those matrices when 7 = 7r/4. Let us define the orthogonal matrix 



O 



[10 \ 

(cos 7 + sin 7) J \[2 (cos 7 — sin 7) / \J~2 
(sin 7 — cos 7) I y/2 (cos 7 + sin 7) A/2 j 



(28) 



Then we easily find that 



O t M r { 1 )0 = M r {ti/A) => O t M v (7) O = M v (tt/4) . 



The matrix 



17' 



/ C12 ^_ S12 

-WV2 c 12 /V2 -1/^2 
V s 12 /V2 c 12 /V2 l/y/2 



where 



diag (l, 



(29) 



(30) 



(31) 



C12 = cos 6*12, and su = sin#i 2 , diagonalizes A^!/(7r/4) 8J. The matrices e ia and e l/3 are 
diagonal phase matrices; the first of them is unphysical, and it can be shown that it 
is of the form given in equation pip. According to equation (|29p. the matrix U which 
diagonalizes M. v (7) is 



U = O U' = e l 



I C12 S12 

— S12COS7 C12COS7 —sin 7 
\ — Si 2 sin7 Ci 2 sin7 cos 7 



(32) 



cf. equation (|23j) . With e 1 ^ = diag (e^ 1 , e 1 ^ 2 , e i/3a ), the physical Majorana phases can be 
chosen to be A = 2 (ft - /3 2 ) and 2 (ft - ft). 

A salient feature of the D4 model — whether softly broken or not — is the fact that 
(Mr) = 0. This translates into jHj 



3=1 



e 2ia sin 7 cos 7 



2 p 2i/h 



s 12 e 



rrii 



+ 



r 2 2i/3 2 
c 12 e 

m 2 



3 2i/3 3 ' 



m 3 



0. 



(33) 



We see that the expression within parentheses is zero irrespective of the value of 7. Hence, 

is, just as in jS], 



the effective mass probed in neutrinoless (3(3 decay, \{m)\ = |(A4„) ee |, 

777.1777.2 



I HI 



7773 



(34) 



Moreover, just as in the original -D4 model, only the normal spectrum m\ < m 2 < 7773 is 
allowed, (mi < 7772 holds by definition, and the solar mass-squared difference is AtUq 



777i 
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For baryogenesis via leptogenesis j.l8j — for reviews, see for instance ^H] — we must 
diagonalize the mass matrix Mr. Denoting the unitary diagonalization matrix by V and 
the heavy- neutrino masses by Mj, one has 

V T M R (7) V = V' T M R (it/A) V = diag (M h M 2 , M 3 ) , (35) 

with V = O V. The matrix which enters the calculation of the CP asymmetry generated 
by the decay of the heavy neutrinos is given by JHj 

R = V T M D M ] D V* = V' T M D M ] D V*. (36) 

The second equality in equation (j3T?j) is justified by the twofold degeneracy of Mo, see 
equation (jUJ). Thus, in the softly broken D4 model only V is relevant for leptogenesis. As 
emphasized before, U' and V coincide with the diagonalization matrices of the original 
-D4 model. Hence, the soft breaking of -D4 does not alter the expression of the lepton 
asymmetry derived in |2()j . and one obtains the same results irrespective of the value of 
7 = #23- The masses of the third-generation neutrinos, both in the light- and heavy- 
neutrino sectors (i.e., both m 3 and M 3 ), do not occur in the expression for leptogenesis— 
which depends only on m 12 , M 12 , #i 2 , and A [20J |2T]. However, m 3 enters leptogenesis 
indirectly, because the Majorana phase A may be expressed as a function of m 3 ; indeed, 
from equation [20], 

m 2 s l2 e + m 1 c 12 = , (37) 

m 3 



hence 

'2 



(mim 2 / 'm 3 ) 2 — c\ 2 m\ — sf 2 m 2 
2c 2 12 s\ 2 m 1 m 2 



cos A = v 1 z/ ^ 2 ^ 1 12 2 . (38) 



Successful leptogenesis in the _D 4 model requires j20j the lightest-neutrino mass m 1 to be 
below 10 -2 eV. In that region equation (}3*Hj) gives a strong restriction on the allowed range 
for mi, since cos A must be larger than —1; using 6\ 2 = 33°, Am^ = 7.1 x 10 _5 eV 2 , and 
Am 2 tm = m\ - m\ = 2 x l(T 3 eV 2 |'3j, one obtains 2.9 x l(r 3 eV < m l < 7.1 x HT 3 eV. 
In order to reproduce the baryon-over-photon ratio i]b, Mi must lie in between 10 11 
and 10 12 GeV, if we assume M 2 3> Mi. The maximum value of rjs is attained for m\ ~ 
4 x 1CT 3 eV, and the experimental value t]b ~ 6.5 x 10 10 |22! can easily be reproduced — see 
figure 1 of 



5 Radiative corrections 

The form of M. v given by equations ©, ©, and holds only at the seesaw scale m R . 
Radiative corrections must be taken into account if one wants to calculate Ai u at the 
electroweak scale. The effective operators relevant in this context are [To"! ITol IT7| 



E E pc*T a K { $c-i{p pL ) t 

a,/3 a,b,c,d 



.cd 



J Vd 



a,/3 



(39) 
(40) 
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where a, b, c, d are SU{2) indices, e = ir 2 is the 2x2 antisymmetric tensor, and the 
are matrices in family space. Note that we may, without loss of generality, enforce the 
condition k^) = Kg^. We want to discuss the renormalization-group (RG) evolution of 
the effective coupling matrices from the seesaw scale ttir down to the electroweak 
scale; the latter may be represented by the Z° mass raz- Our aim is to estimate a possible 
modification of equation ([2fip by the RG evolution. We denote the dimensionless variable 
of the RG by t; the RG evolution goes from to = ^{mn/mz) to t\ = 0. The initial 
conditions for the RG equations are 



K W) (*„) = for all other (ij). 



(41) 



In equation (j4~Tj) we have used the fact that in our model only the Higgs doublet <pi has 
Yukawa couplings to the right-handed neutrinos, and gives thereby rise to Mr>. 

The symmetry Tj 2 aux \ which is preserved throughout the RG evolution since it is 
only broken spontaneously at the electroweak scale, implies that there is a symmetry 
(pi — > —0i in the evolution of the operators in equation flJOJ). Since at the high scale only 
On is present, which is invariant under Tj^^, it follows that the operators Ou and O13, 
which are odd under zjf , remain zero at all scales. At the electroweak scale we will 
then have 

M v = £ (h) + v 2 v 3 [k^ ( tl ) + (h)] . (42) 

The matrix M. v at the electroweak scale will not in general yield U e s = 0. Indeed, 
using U T M. V U = diag (m 1 , m 2 , m 3 ), we find that U e3 = if and only if there is a vector u 
with a zero first entry such that Ai u u = zu*, where z is a complex number and m 3 = \z\. 
This means that the symmetric matrix A4 U satisfies 

(M,)'-I(M,)J 



{M V )^ T = {M v )* e , (M u ) eT {M v ) m - (MJ^ (M„T eT {M V ) TT . 

(43) 

This condition will in general not be satisfied by the matrix M. u at the electroweak scale 
given in equation ()42|1. but it is satisfied by the matrix M. v at the seesaw scale given by 
equations ®, ©, and Q. 

Denoting the set of all matrices k^' by k, the differential equation for their RG 
evolution can symbolically be written as 



2 d «W 



1671 dt =Lt [K{t)] ' (44) 

where L t is a linear operator acting on the In L tl the Yukawa and gauge coupling 

constants appear in second order and the quartic Higgs couplings in first order. The 
solution of equation (|44|) is formally given by the series 



n(t) = K(t ) + ~^f dt'L t , [«(*„)] + (-M f dt' f dt"L t , [L fl [K(t )]\ + ■ • ■ (45) 

1 7T J to \ 1 7T / J to J to 
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Figure 1: A typical vertex correction which introduces flavour dependence in the renor- 
malization of the operators Oij of equation (|4()jl . The relevant Yukawa-coupling matrices 
are indicated. 

Flavour dependence in k^' (the dot denotes the derivative relative to t) can only be 
introduced by the Yukawa couplings. One can easily convince oneself that only vertex 
corrections of the type depicted in figure [T] need a closer look at. 2 The Yukawa-coupling 
matrices of the Higgs doublets (f)j to the right-handed charged leptons are given by the 
second line of the Yukawa Lagrangian in equation (J2J); they are 

Y 1 = diag (</ 3 *, 0,0), 

Y 2 = diag (0,yl,yt), (46) 
Y 3 = diag (0,yl-y* 5 ). 

One finds from figure[T]that k^ + k^ obtains a contribution YJ fc YjYk + k^Y^Y^ 
plus the transposed term. 

A non-zero is induced from K^(to) in two steps. Firstly, vertex corrections from 
the quartic terms 

X' (4<p 2 y + X" (4^ (47) 
in the Higgs potential induce non-zero matrices and via 

l 6vr 2^ A / K (U) (48 ) 

dt y ' 

IQ^^-l „ AV 11 ). (49) 
dt 



2 All the types of graphs relevant for the computation of the full RG equations for the «W) are depicted 
in [E3E7]. 
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Secondly, from a vertex correction of the type in figure ^ one obtains 



16**1 

dt 



+ K< 32 > 



K (22)y 3 ty 2 + Y?Y*K<& + /€ {33 V 2 V 3 + Y 3 T Y*K^\ (50) 



Equations (jlEJ), P^jl - and (f5T)|) tell us that, in the expansion of equation (jlKjl . the 
matrices ft*- 23 -* and k^ 32 ^ appear first at second order. Since 

Y*Y 2 = MY 3 V = y*y B diag(0, +1,-1), (51) 



the vector w 3 of equation (|23|) will not be eigenvector of z^ 23 ) + k^ 32 \ One can check 
that, at second order in the expansion of equation (|45p. u 3 is still an eigenvector of the 
K {n) ^ _ ^ 2, 3). Therefore, we estimate that the terms in A4 U which are responsible for 
U e3 7^ will typically be suppressed by a factor of the order of 

1 /, m R \ 




In^l 1^1 | A '."|. (52) 



(167T 2 ) 2 V rrizJ 

In this equation we have taken into account that there are two integrations over the 
interval of length t Q and that the relevant Yukawa couplings are of the order m T /\v2\- 
Numerically, using t ~ 10 and even if we take the quartic Higgs couplings to be of order 
one, the suppression factor is at least 10 -4 . Therefore, to a good approximation we will 
still have £7" e3 = at the electroweak scale. 3 



6 Conclusions 

In this paper we have modified the D4 model of [S] by allowing a soft breaking of the 
subgroup of D4, i.e. by breaking the /i-r exchange symmetry softly — see equation ((TJ). 
An essential ingredient of the -D4 model is the -D4 doublet (X11X2) of real scalar gauge 
singlets, with VEVs and masses of the seesaw scale rriR. It is the VEVs of those scalars 
which induce lepton mixing. In particular, one has 

tan# 23 = -^. (53) 
(Xi)o 

(tr) 

Requiring the soft breaking of Z 2 to occur exclusively through terms of dimension two 
in the Lagrangian, we find that the term of equation (fTSj) appears in the scalar potential. 
Whereas without soft breaking the VEVs of xi an d Xi are (almost) equal, the soft breaking 
of 7L 2 induces a deviation from this equality. Thus, the strength of the soft-breaking 
term in equation (|18j) determines the deviation of #23 from 45°, small deviations being 
natural in a technical sense. 

We have demonstrated that the soft D4 breaking has no effect on the neutrino mass 
spectrum, on the effective mass | (m) | in neutrinoless (3(3 decay, and on leptogenesis. Just as 
in the original -D4 model, we predict a normal spectrum m x < m 2 < m 3 , < 10~ 2 eV, 

and we find that leptogenesis is successful when mi ~4x 10~ 3 eV. 

3 The details of the RGE for multiple-Higgs-doublet models will be published elsewhere. 
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Like in all multi-Higgs models where the Yukawa-coupling matrices are diagonal and 
lepton mixing originates solely in Mr, neutral flavour-changing interactions are suppressed 
|13j : on the one hand, the amplitudes for e.g. yT — ► e~7 and Z Q — > e + fi~ are propor- 
tional to m^ 2 , hence highly suppressed; on the other hand, a decay like \i~ — > e~e + e~ 
is suppressed by Yukawa couplings but might fall in the discovery range of forthcoming 
experiments. 

The -D4 model — whether softly broken or not — predicts # 13 = 0. We have estimated, 
through renormalization-group methods, the radiative corrections to this relation to yield 
#13 ~ 10 -4 or smaller — the same estimate also applies to the model of jHj. The solar 
mixing angle is not predicted in the D4 model; without finetuning it will be large. 

In summary, we have constructed an extension of the Standard Model based on the 
non-Abelian family symmetry group D4. Our model features a vanishing #13 together 
with non-maximal atmospheric mixing. The model is in agreement with all the existing 
data on neutrinos and allows successful leptogenesis. 
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